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Abstract 

We study the influence of the softness of the interparticle interactions on the fragility of a glass 
former, by considering three model binary mixture glass formers. The interaction potential between 
particles is a modified Lennard- Jones type potential, with the repulsive part of the potential varying 
with an inverse power q of the interparticle distance, and the attractive part varying with an inverse 
power p. We consider the combinations (12,11) (model I), (12,6) (model II) and (8,5) (model III) 
for (q,p) such that the interaction potential becomes softer from model I to III. We evaluate the 
kinetic fragilities from the temperature variation of diffusion coefficients and relaxation times, 
and a thermodynamic fragility from the temperature variation of the configuration entropy. We 
find that the kinetic fragility increases with increasing softness of the potential, consistent with 
previous results for these model systems, but at variance with the thermodynamic fragility, which 
decreases with increasing softness of the interactions, as well as expectations from earlier results. 
We rationalize our results by considering the full form of the Adam-Gibbs relation, which requires, 
in addition to the temperature dependence of the configuration entropy, knowledge of the high 
temperature activation energies ino rder to determine fragility. We show that consideration of 
the scaling of the high temperature activation energy with the liquid density, analyzed in recent 
studies, provides a partial rationalization of the observed behavior. 



I. INTRODUCTION 



The temperature variation of relaxation times, viscosity and diffusion coefficient in glass 
forming liquids upon approaching the glass transition has been studied for a wide variety 
of substances. Near the glass transition, these quantities show a rapid increase, but with 
a rate of change that is different for different substances. The rapidity of rise of relaxation 
times near the glass transition has been quantified by "fragility", introduced and analyzed 
extensively by Angell which has proved to be useful in organizing and understanding the 
diversity of behavior seen in glass formers. Fragility has been defined in a variety of ways. 
Two of the popular definitions are in terms of the "steepness index" m, and the fragility 
defined using Vogel-Fulcher-Tammann (VFT) fits to viscosity and relaxation time data. 

The steepness index of fragility is defined from the so-called Angell plot as the slope (m) 
of logarithm of the viscosity (77) or relaxation time (r) at T = T g , with respect to the scaled 
inverse temperature T g /T where T g is the laboratory glass transition temperature: 

m = ( d -^T) (1) 

We refer to the fragilities defined from transport quantities and relaxation times as kinetic 
fragilities, to be distinguished from thermodynamic fragilities defined later. A kinetic fragility 
may also be defined from a VFT fit of the relaxation times, 



r(T) = r exp 



f 



(2) 



which defines the kinetic fragility Kypr and the divergence temperature Typr- 



Despite considerable research effort 



lMl7j , and the observation of many empirical corre- 
lations between fragility and other material properties, a fully satisfactory understanding of 
fragility hasn't yet been reached. Such understanding has been sought, broadly, along two 
lines. The first is a conceptual understanding of fundamental quantities that may govern 
fragility. An example of this kind is the use of the potential energy landscape approach in 
combination with Adam Gibbs (AG) relation [18] between relaxation time and configuration 
entropy [Eq. |3] to relate features of the energy landscape of a glass former to the fragility. 
The Adam-Gibbs relation 

r(T) = r exp(^|£) (3) 



relates the temperature dependence of the relaxation times to the temperature change in the 
configuration entropy S c , where 5{i is an activation free energy for particle rearrangements, 
and S* is the configurational entropy of cooperatively rearranging regions invoked in Adam- 
Gibbs theory If A = dfiS*^ 1 has no significant role to play in determining the fragility 
of a substance, it is the temperature variation of TS C that dictates the fragility. If the 
T-dependence of S c is given by 

TS e = K T (I- - l) , (4) 

the Adam-Gibbs relation yields the VFT relation, with the identification Kyp? = Kt/A, 
Tvft — Tk- Thus, Kt is a thermodynamic index of fragility. 

In what follows, we use Eq.s [2] and [3] which describe our simulation data well, as we 
demonstrate. However, our discussion does not depend crucially on the strict validity of the 
VFT temperature dependence near the glass transition, or the divergence of relaxation times 
at finite temperature; both these features have been questioned by various investigations and 



alternative forms to the VFT temperature dependence have been proposed 



In potential energy landscape approach 22|, [231] configuration entropy is associated with 
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the number of local potential energy minima or inherent structures (IS) [24(, and can be 
computed in terms of parameters describing the energy landscape [3J]. Hence thermodynamic 
fragility can be understood in terms of parameters of the potential energy landscape, namely 
the distribution of inherent structures and the dependence of the vibrational or basin entropy 
corresponding to inherent structures on theie energies. Although the exact temperature 
dependence of the configuration entropy depends on detailed properties of the distribution 
of inherent structures, and is not a constant even in the simplest case, such analysis 
does yield insight into the relationship between the energy landscape features and fragility. 
To a first approximation, the broader the distribution of energies of inherent structures, 



the larger the fragility of a glass former [3]. Going beyond such analysis, one needs to 



also understand the 



aehavior of the prefactor A, which is related to the high temperature 



activation energy 



251 ] . To the extent that the Adam-Gibbs relation quantitatively 
describes the temperature dependence of the relaxation times, such analysis provides a 
route to a fundamental understanding of fragility in terms of the phase space properties of 
a substance. 

However, such a conceptual understanding does not directly address the dependence of 



fragility on specific, controllable material properties, an understanding that is desirable from 
the perspective, e. g., of materials design. The investigation of the dependence of fragility 
on the nature of molecular architecture and intermolecular interactions defines therefore 
a second distinct line of investigation, which has been pursued by various groups. For 
example, Dudowicz, Freed and Douglas 



10 



111 ] have investigated the role of backbone and 



side group stiffness in determining the fragility of polymer glass formers. In another recent 
example, from an experimental investigation on deformable colloidal suspensions, Mattsson 



et al 



16 



171 ] suggested that increasing the softness of the colloidal particles should decrease 



the fragility of the colloidal suspensions, and that such a principle should be more generally 
applicable. Indeed, this conclusion is consistent with that of Douglas and co-workers 10] 
that the ability to better pack molecules leads to lower fragilities. In energy landscape terms, 
one may understand this conclusion as implying that molecules that pack well together will 
have narrower distributions of inherent structure energies. 
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FIG. 1. Comparison of interaction potential V a p without truncation for the three different poten- 
tials used in the present study, r™™ are the positions of the minima of the interaction potentials. 



The influence of the softness of interaction on the fragility was also investigated some time 
ago via computer simulations of model glass formers by Bordat et al js, 9]. They considered 
a binary mixture of particles interacting via generalized Lennard Jones potentials, of the 
form 

nr) = (? ± ri( p( 7 ), -4>') < 5 » 

for combinations of the exponents (q,p) of repulsive and attractive parts of the potential 
(12,11), (12,6) and (8,5). These combinations, corresponding to models labeled I, II and III, 
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FIG. 2. Comparison of the interaction potential V a p with truncation for the three different poten- 
tials, r™™ are the positions of the minima of the interaction potentials. 



have decreasing curvatures at the minimum of the potential, and thus increasing softness. By 
evaluating the kinetic fragility of these models (the steepness index defined above), Bordat 
et al found that increasing softness of the interaction potential increases the kinetic fragility 



a. 



The trend found by Bordat et al therefore is apparently not consistent with expectations 
arising from the other studies mentioned, although the nature of the changes in the inter- 
actions considered are not strictly the same. In order to understand better the relationship 
between the nature of the intermolecular interactions and fragility, in the present work we 
calculate the kinetic fragility Kvft using computer simulation data of the diffusion coef- 
ficient, and relaxation times obtained by a number of different means. We also calculate, 
using the procedure in 0, Q, the configuration entropy, from which we calculate a 
thermodynamic fragility {Kt). We find that these two fragilities show opposite trends, with 
the kinetic fragility increasing with softness, and the thermodynamic fragility decreasing 
with softness. In order to understand this apparent disagreement, we must consider the 
full form of the Adam-Gibbs relation, including terms that relate to the high temperature 
activation energy. We present our analysis along these lines below. We focus our analysis 
here on the role of a specific feature of the interaction potential, namely the softness, for 
reasons stated above. However, fragility in principle depends on a number of parameters 
that describe a glass former, which may include pressure, density etc.. While our study 
implicitly includes those factors that are affected by a change in the interaction potential 



(keeping other parameters fixed), we do not attempt here a comprehensive analysis of all 
factors that may influence the fragility of a glass former. Related questions concerning the 
change in structure, dynamics and thermodynamics in a glass forming liquid upon tuning 
the interaction potential have been addressed in 



The paper is organized as follows: In Section II we summarize the computer simulation 
details. In Section III we describe the methods used for evaluating the various quantities of 
interest. In Section IV we present our results and a discussion of the results, and Section V 
contains our conclusions. 



II. SIMULATION DETAILS 

We have studied a 80:20 binary mixture of modified Lennard Jones particles in three 
dimensions. The interaction potential is of the form given above in Eq. [5] with a truncation 
that makes both the potential and force go to zero smoothly at a cutoff distance r c . The 
potential with the truncation is given by 



\r i \ e ° l l 3 
Vap[r) = — ^ 

q-p 



„mm „mm 

p( ^L )9 _ g( ^L ; 



+ci a pr + c 2a p,r < r ca)3 
= 0, otherwise (6) 

where a, f3 G {A,B}. r™™ = 2^a al 3 and e a p are respectively the position and the value of 
the minimum of the pair potential. The correction terms Ci a p, C2 a p are determined from the 
conditions : 



V a p{r ca p) = 
dr 



= (7) 

fee /3 

The energy and size parameters and a a/ 3 correspond to those of the Kob-Andersen 



binary Lennard- Jones model 28]. Units of length, energy and time scales are oaa^aa 
and W ^g^ - respectively. In this unit, €ab = 1-5, €bb = 0.5, <jab = 0.80, obb = 0.88. 
The interaction potential was cutoff at 2.5er a| g, The three different models (12,11), (12,6) 
and (8,5) are shown without and with cutoff in Fig. s CD and El Molecular dynamics (MD) 
simulations were done in a cubic box with periodic boundary conditions in the constant 
nummber, volume and temperature (NVT) ensemble. The integration time step was in the 



range dt = 0.001 — 0.005. Temperatures were kept constant using an algorithm due to Brown 
and Clarke u 



Simulations were done in the temperature range T G [0.85,5] for (12, 11); 
T G [0.45, 5] for (12, 6) and T G [0.23, 5] for (8, 5) model respectively. System size were 
N = 1500, Na = 1200 (N = total number of particles, Na = number of particles of species 
A) and the number density was p = 1.2 ([8j, see also Fig. [3]). For all models, one sample 
per state point above the onset temperature (described below) and three to five samples 
per state points below the onset temperature were used with runlengths > IOOtq, (r a is the 
relaxation time, described below). 




FIG. 3. Pressure vs. density for inherent structures (IS). The density minimum for IS pressure 
occurs at p = 1.04, 1.09, 1.18 respectively for models I (12, 11), II (12, 6), and III (8, 5). This density 
defines the lower bound for simulations of the system in the homogeneous liquid state. 



III. METHODS 

In this section, we describe the various quantities that have been calculated and the 
methods employed for such calculations. 



A. The a relaxation time 

The following measures have been used to extract a relaxation times: 
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1. Diffusion coefficient (Da) from the mean squared displacement (MSD) of the A type 
particles. 

2. Relaxation times obtained from the decay of overlap function q(t) using the definition 
q(t = r a ,T)/N = 1 /e. The overlap function is a two-point time correlation function 



of local density [3CH34j which has been used in many recent studies of slow relaxation, 



and is defined as: 



< q(t) > = < J dfp(r, t )p(f, t + t )> 



N N 
i=l j=l 

Here the averaging over time origins to is implied. The overlap function naturally 
separates into "self" and "distinct" terms: 

N 

< q(t) > = <J2 fifth) ~ r t (t + t )) > 

+ <J2J2 6 ft^ to ) - r 3 (t + to)) > 

In our work, we consider only the self part of the total overlap function (i.e. neglect 



the i 7^ j terms in the double summation), based on the observation 32] that the 
results obtained from the self part are not significantly different from those obtained 
by considering the collective overlap function. Thus we use 



< q(t) > w < Sftih) - n(t + t )) > 

i=l 

Further, for numerical computation, the 8 function is approximated by a window 
function w(x) which defines the condition of "overlap" between two particle positions 
separated by a time interval t: 

N 

< q(t) >^< J2^(ft(to)-ri(to+t)\) > 
i=i 

w(x) = l,x < a implying "overlap" 
= otherwise 



S 



The time dependent overlap function thus depends on the choice of the cutoff parame- 
ter a, which we choose to be 0.3. This parameter is chosen such that particle positions 
separated due to small amplitude vibrational motion are treated as the same, or that 
a 2 is comparable to the value of the MSD in the plateau between the ballistic and 
diffusive regimes. 

3. We have also studied the "susceptibility" x±{t)-> defined in terms of the fluctuations in 
the overlap function as 

X*(*) = ^ {(lit) 2 ) - (lit)) 2 ) (9) 
This quantity can be written as an integral of a higher order, four point correlation 



function g±{r, t) 



30 



3ll . |33| widely studied in the context of dynamical heterogeneity: 

g 4 (f,t) = (p(0,0)p(r,0)p(0,t)p(r,t)> - 
(p(0,0)p(f,0))(p(0,t)p(f,t)) 
XA (t) = [ drg A {r,t) (10) 



The characteristic time ta(T) at which the fluctuation (Xiit)) is maximum is taken as 
a measure of relaxation time. 

4. Relaxation times obtained from the decay of the self intermediate scattering function 
FJk.t) using the definition FJk.t = T a ,T) = lie at k ~ The self intermediate 

Train 

scattering function is calculated from the simulated trajectory as: 

1 N 

F s (k,t) = ^<X>p (-«£• ft® -n(o)))) (ii) 

Since the relaxation times from q(t), Xiif) ari( l Fs{k,t) behave very similarly, we discuss 
further only the time scale obatined from q(t). 



B. Characteristic temperature scales 

Dynamics of fragile glass forming liquids show characteristic cross-over from high tem- 
perature Arrhenius behaviour to low temperature super Arrhenius behaviour at some char- 
acteristic temperature. At this temperature, systems also show cross-over from a landscape 
independent high temperature regime to a landscape influenced low temperature regime 

9 



TABLE I. Characteristic Temperatures 



Quantity (12,11) (12,6) (8,5) 



- onset 



1.27 0.9 0.42 



T c from {^y 1 0.77 0.42 0.22 

T c from q(t) 0.77 0.42 0.23 

T VFT from (^r 1 0.59 0.32 0.17 

T VFT from q(t) 0.55 0.29 0.16 

T K 0.54 0.28 0.16 



35 



361 ] . We denote this temperature as onset temperature T onset . We report the estimates 
from inherent structure energies in Table [I] As the temperature is further lowered, mode 
coupling theory predicts divergence of relaxation time r as t(T) ~ (T — T c ) -7 which defines 
the mode coupling divergence temperature T c which we estimate from both relaxation time 



and diffusion coefficient (in the form (Da/T) -1 ) 28]. Similarly relaxation times apparently 
diverge at a second characteristic temperature which we estimate from VFT fits and denote 
as TypT- Further, configuration entropy becomes zero on extrapolation at a characteristic 
temperature (Eq. Hj) known as Kauzmann temperature (Tk). The AG relation (Eq. [3]) 
predicts that these two temperatures (T VFT and T K ) to be same. Although we use func- 
tional forms that have a temperature of vanishing S c and diverging relaxation times, these 
are employed as useful descriptions of the data, without any implied assertion of the ex- 
pected behavior at temperatures lower than the ones we study. The The values of different 
characteristic temperatures for different potentials are tabulated in Table [D 



C. Configuration entropy 

Configuration entropy (S c ) per particle, the measure of the number of distinct local energy 
minima, is calculated 26[ by subtracting from the total entropy of the system the vibrational 
component: 

S C (T) = S total (T) - S vib (T) (12) 

The total entropy of the liquid is obtained via thermodynamic integration from the ideal 
gas limit. Vibrational entropy is calculated by making a harmonic approximation to the 

10 



- O Model I (12,11) 
6 - □ Model II (12,6) 
O Model III (8,5) 




FIG. 4. Temperature dependence of TS C for the studied models to determine the Kauzmann 
temperature. T K = 0.54,0.28,0.16 respectively for models I (12, 11), II (12,6), and III (8,5). The 
value of Tk from the extrapolated crossing of bulk and basin entropies vs. temperature reported 



m 



270 is T K = 0.2976 and in 



371 ] is Tk ~ 0.29. The Tk values obtained from this plot is used to 



determine the thermodynamic fragility in Fig. [5j 




T/T, 



FIG. 5. Determination of the thermodynamic fragility from the relation TS C = Kt(jt- — 1) for 
the studied models. Kt is the slope of the linear fit shown. Tk is the temperature at which 
S c = 0, obtained from the linear fit shown in Fig. [4j Thermodynamic fragility (Kt) values are 
0.551,0.323,0.211 for models I (12,11), II (12,6), and III (8,5) respectively 
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potential energy about a given local minimum [22l . l23l l26l . |27| |. The procedure used for 
generating local energy minima, and calculating the vibrational entropy is as outlined in 

We have also computed the configuration entropy density S c (ejs) = k B \iafl(eis) where 
Q(ejs) is the number density of inherent structures with energy eis and to a good ap- 
proximation may be described by a Gaussian. Equivalently, S c (eis) can be described by a 
parabola 



S c (e IS ) = a- ™ ; is) (13) 

The parameter a denotes the peak value of S c (eis) which occurs at energy e° Is . S c (eis) is 
zero at ejs = e% i °\f^- Thus a^fa is a measure of the spread of S c {ejs). We denote the 
lower root e Q IS — o^fa by e^ n . 

In the harmonic approximation to vibrational entropy, the average value of IS energy 
sampled by a system at a given temperature < e/s > (T) is predicted to be linear in inverse 
temperature (3 = 1/T: 

(eis)(T) = < ej s)(oo) - ^ (14) 

where (e/s)(oo) is the extrapolated limiting value of (ejs) at high temperatures. These 
parameters which characterizes the potential energy landscape are tabulated in Table [III for 
different potentials. 



IV. RESULTS AND DISCUSSION 

The results from the simulations concerning the thermodynamic and kinetic fragility 
estimates are presented below. 



A. Thermodynamic fragility 

As described in the Introduction, we may define a thermodynamic fragility as the slope 
of TS C (T) vs. T/T K . Fig. 0] shows that indeed, TS C (T) varies linearly with temperature, 
which allows us to define Tk- The Tk values for the different potentials are listed in Table 
U Various quantities related to the distribution of inherent structure energies are listed in 
Table II for later use. Thermodynamic fragility as defined in Eq. H] is computed from 

12 



the slope of TS C vs. T/Tk is found to decrease as the softness of the interaction potential 
increases, as shown in Fig. [51 Such behavior is in line with expectations, e. g. from lOl. Il6|. 



B. Kinetic fragility 

Kinetic fragility (Kvft) is estimated by fitting to the VFT form Eq. [2] the diffusion 
coefficients and relaxation times. 

In Fig [6] (top panels), we show the Arrhenius plot of the diffusion coefficients and 
relaxation times from q(t), plotted against Tk/T. The VFT divergence temperatures Tvft, 
obtained from VFT fits to the data for temperatures below the onset temperature, are found 
to be close to and are listed in Table I. The middle panels of Fig [HI show Arrhenius 
fits to high temperature data (above the onset temperature), from which activation energies 
Eq (such that t(T) = tq exp(£'o//csT) are obtained. These are listed in Table III, and will 
be discussed later. In the bottom panels of Fig [61 we show the effective activation energy 
defined as E{T) = A; B Tln(r(T)/r ) scaled by E (similarly for D A /T), plotted against 
k B T/E . 

We note in the passing that for model (12, 6) the proportionality Eq ~ 6T C [10j is rea- 
sonably well satisfied. However, the ratio E /T c decreases from ~ 7 to ~ 5 as softness 
increases. 

Next, we calculate the kinetic fragilities Kvft, from diffusion coefficients and relaxation 
times, using the divergence temperature T VFT obtained with T VFT as a fit parameter, as 
well as using Tr estimates from the configuration entropy as the divergence temperatures. 
The corresponding kinetic fragilities, labeled K VFT and K\} FT , are listed in Table IV, along 
with the thermodynamic fragilities K^. We find that the kinetic fragilities increase as the 
softness of the interaction potential increases, thus showing a trend that is opposite to that 
of the thermodynamic fragility. 

C. Adam Gibbs relation and fragility 

In order to understand this discrepancy, we consider again the Adam-Gibbs relation, 
which relates the kinetic and thermodynamic fragilities. Comparing Eq. [2J, Eq. [3] and Eq. 
HI we note that the relationship between the kinetic and thermodynamic fragilities that we 
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may deduce assuming the validity of the VFT and the AG relations is 

K VFT = K T /A (15) 

and we expect at least the same trend in the two fragilities under the assumption that the 
term A does not substantially alter the proportionality between kinetic and thermodynamic 
fragilities. To assess the degree to which this is true in our models, we show in Fig. [7] the 
Adam-Gibbs plots of the diffusion coefficient and relaxation times. These plots show that 
the coefficient A, obtained from the slopes (and listed in Table III), indeed varies from one 
model to the other, decreasing as the softness increases. Thus, the ratio Kt/A shows the 
opposite trend, increasing as the softness increases. 

We next attempt to understand the dependence of the Adam-Gibbs coefficient A on the 
softness of the interaction. First we consider the high temperature Arrhenius behavior of 
relaxation times, in terms of the Adam-Gibbs relation. Such Arrhenius behavior can be 
expected if the configuration entropy effectively becomes a constant, in which case, the high 
temperature activation energy will be given by 

E = Ak B /S CoD (16) 

However, the asymptotic high temperature configuration entropy is difficult to assess di- 
rectly, as the various available approaches to computing the basin entropy do not work well 
in this regime(see e. g. 36]). We thus use the following procedure: First, we determine 
directly from simulations the high temperature limit of the inherent structure energies, e 1 }™ 
(see Fig. Then, we use the extrapolation of the dependence of the configuration entropy 
S c on the inherent structure energy e/s obtained below the onset temperature to obtain the 
high temperature limit of the configuration entropy, S^e}™), which do not vary appreciably 
with softness of interaction, and are listed in Table II. Table II also lists S c (oo), the infinite 
temperature value of S c obtained by extrapolating Eq. H] to infinite temperature, a proce- 
dure that is not justified at temperatures above the onset temperature. Using these S c {e 1 }^) 
values, and the activation energies Eq shown in Table III, we obtain estimates for the AG 
coefficient 

A est = E S c (e l i™)/k B (17) 

which are shown in Table III. We note in Table III that Eq values decrease strongly as 
the softness of the interactions increases, and with a corresponding moderate increase of 

14 



Sd&j'g'), our estimates of A est agree rather well with the values obtained directly from the 
Adam-Gibbs plots. We now designate the thermodynamic fragility estimates obtained by 
considering the full form of the Adam-Gibbs relation as Kag — Kt/A, and list them along 
side the thermodynamic and kinetic fragility estimates in Table IV. As expected from the 
above discussion, the "Adam-Gibbs" fragility estimates (K AG in Table IV) agree rather well 
with the kinetic fragilities. 

Although the above picture provides a consistent description of the fragilities from kinetic 
and thermodynamic data, a question remains regarding the variation of the high temperature 
activation energy Eq with the softness of the interaction potential. To seek some insight into 
this question, we consider work in recent years concerning the scaling of the temperature 
dependence of dynamic and thermodynamic quantities at different densities j^, Q, 38 1 . It 
has been shown by many groups that a scaled variable p 7 /T, where p is the density, captures 
the density variation of properties in many liquids. The exponent 7 can easily be shown 
to be n/3 for inverse power law potentials, where n is the power of the inverse power law, 
but even for other liquids, an effective 7 has been shown to be derivable by considering the 
correlated fluctuations of potential energy and the virial [25]. The exponent 7 is obtainable 
as the ratio of fluctuations. Although such a ratio is state point dependent, a "best fit" value, 
typically obtained from high temperature state points, has been shown to effectively describe 
the scaling of properties at different densities. Since we do not perform a full analysis of the 
density dependence here, we do not estimate the best value of 7 but instead use the value at 
twice the onset temperature as an indicative value. Fig. [9] shows the fluctuation data from 
which the 7 value is obtained, and the temperature variation of the exponents. The values 
of 7 we use are shown in Table II. 

Based on the above considerations, we should expect the high temperature activation 
energies to be proportional to p 1 . Accordingly, we obtain estimates of the activation energy 
in the form Eq = Eqqp 1 . These values, shown in Table III, have a weaker temperature 
dependence than the directly evaluated E , and correspondingly, the fragility estimates 
obtained (shown in Table IV), while showing a smaller decrease with softness, nevertheless 
decrease with increasing softness of interaction. A further analysis is needed, therefore, 
to elucidate the relevance of these considerations to evaluating the variation of the high 
temperature activation energy. 
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V. CONCLUSIONS 



We have studied the effect of the softness of the interaction potential on fragility in three 
model glass formers. We find that the kinetic fragility obtained from diffusion coefficients 
and relaxation times increases with increa sing softness of the interaction potential, contrary 
to expectations based on earlier studies [10l . Il6j . On the other hand, a thermodynamic 
fragility obtained from the temperature variation of the configuration entropy decreases 
with increasing softness of the interaction potential. By taking into consideration the model 
dependence of the high temperature activation energy, in addition to the temperature de- 
pendence of the configuration entropy, we define an "Adam-Gibbs" fragility whose model 
dependence accurately captures the variation of the kinetic fragilities that we find. An at- 
tempt to rationalize the model dependence of the high temperature in terms of the scaling 
of properties with respect to density is encouraging but fails to fully explain the observed 
decrease of the fragility with increasing softness of the interaction potential. 
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FIG. 6. Top row: Inverse diffusion coefficient and relaxation time from overlap function vs. scaled 
inverse temperature ^F. Lines through the data show VFT fits to the data below the onset 
temperature. T% estimated from fig. H] are used as the divergence temperatures in the VFT fits. 
Middle row: Arrhenius fits to high temperature data of inverse diffusion coefficient and relaxation 
time from overlap function to determine high temperature activation energies E$. Bottom row: 
effective activation energy E(T) (see text) scaled by Eq, plotted against ksT/EQ. 
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FIG. 7. Adam Gibbs plots for the inverse diffusion coefficient of A particles and relaxation time 
from overlap function, for the three models studied. The activation energy parameter A in Eq. [3l 
obtained from the slopes of the data shown, is tabulated in table IIII1 
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FIG. 8. Temperature dependence of the inherent structure energy eis for the studied models 
shifted by the corresponding high temperature limiting values ejo 1 for clarity. The values of ef™ 
are -6.003,-6.886,-7.191 for models I (12,11), II (12,6) and III (8,5) respectively. 
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FIG. 9. Determination of the density-temperature scaling exponent 7 from the correlation between 
instantaneous potential energy (U) and virial (W). 71 = -^jpp 72 = ^f^) 2 ) ' 73 = $£wau) 
where AU = U — (U) and AW = W — (W) represent fluctuations about mean of potential 
energy and virial respectively. The left panel shows the correlation between energy and virial at 
temperatures « 2T onset , with straight line fits (AW) = 72(2T onset )(AC7). The right panel shows 
the temperature dependent values of 7 for the studied models. 
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